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(57) Abstract: The present invention concerns a method and an arrangement for reducing noise in a substance raw data image 
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containing noise. The present invention also relates to a method and an arrangement for determining a measure of at least one 
substance. Each image is associated to one detection microcavity and one substance that is present in the detection microcavity of 
""--^ a microfluidic structure of a microfluidic device. The substance raw data information is obtained by use of a detector arrangement 
capable of scanning a search area associated to a detection area, which is associated to the detection microcavity. The method may 
comprise following steps in any order for removing different noise signals contributions step-by-step in any order in a substance raw 
data image : - reducing background radiation (step a); - reducing peak disturbance (step p); - locating the detection area within the 
search area (step %); - moving binary artefacts (step A); - removing unwanted areas of the detection area (step e); - applying default 
1^" detection area in noisy images (step ())). 
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Method and arrangement for reducing noise in a 

microfluid device 

5 BACKGROUND OF THE INVENTION 

I. Field of the Invention 

The present invention generally concerns the field of microfluidic devices. More 

particularly, the present invention concerns a method and an arrangement for reducing noise in a 
substance raw data image containing noise. The present invention also relates to a method and an 
10 arrangement for determining a measure of at least one substance. The substance raw data image 
has been obtained from a search area that comprises a detection area associated with a detection 
microcavity that is part of a microchannel structure of a microfluidic device detectable signals. 

The substance raw data image reflects from a substance that is present in the detection 
microcavity after one or more liquid aliquots have been processed through the microchannel 
15 structure that comprises the detection microcavity. See further below. 

II. Related Art 

All patents and patent applications cited in this specification are incorporated by reference in 
their entirety. 



20 A. Background Publications 

Detector arrangements for measuring radiation signals from individual detection 

areas on a circular substrate have been described in a number of previous publications. See for 
instance: EP 392475 (Idemitsu Petrochemical Co, Yamaji Kazutaka et al) US 5994150 (Imation 
Corp, Challener et al\ US 5892577 (The University Court of the University of Glasgow. 
25 Gordon), WO 9721090 (Gamera, Mian et al\ Duffy et aL "Microfabricated Centrifugal 
Microfluidic systems: Characterization and multiple Enzymatic Assays" (Anal. Chem. 71 (1999) 
4669-4678), WO 0040857 (Amersham Pharmacia Biotech AB, BjOrkesten et al) 

B. Background Technology and Problems 

The present invention belongs to the field of miniaturization of processes which 

30 comprise sample treatment, assay protocols, chemical and/or biochemical synthesis etc., within 
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medicine, chemistry, biochemistry, molecular biology and the like. At present one important 
goal within this field is to reduce the costs for these processes, for instance to reduce the amount 
of reagents needed per assay, reduce time per assay, etc. One route has been to increase the 
degree of parallelity, for instance by integrating as many as possible of similar process runs in 
5 one and the same device in order to carry out all the runs in parallel. At present, large numbers of 
research groups and companies are involved in developing technology that will solve the 
numerous problems encountered. 

One problem is related to the optimal way of configuring the detector in relation to 
the microdevice used for performing the processes while maintaining an acceptable sensitivity 
10 and reproducibility. This problem may become particularly pronounced if the measuring step is 
performed by continuously moving the detector unit and the detection areas of a micro fluidic 
device relative to each other during the measurement operation. 

Another problem arises if the microfluidic device is in the form of a disc which is 
skewed because then it becomes difficult to maintain the optical focus in the right position 
15 relative to the detection areas. Without proper arrangement skewed discs will reduce sensitivity. 
This problem in particular applies to discs made of plastic material. 

Another problem is related to maintaining an acceptable sensitivity and 
reproducibility when changing sample volumes from the pi-range to the nl-range and performing 
the process protocols with a high degree of parallelity within the same device. The inventors 
20 have found that under these circumstances the materials from which the microfluidic devices are 
fabricated and the various treatments during the manufacturing and conditioning of the devices 
easily introduce signal artifacts that are of the same kind and of comparable or larger size as the 
desired signals. 

During recent years it has become popular to fabricate microfluidic devices in 
25 plastic material. This kind of material is typically highly fluorescent ("auto-fluorescent") with 
emission wavelengths covering most of the wavelengths normally utilized in fluorescent 
measurements. Compared to microtitre wells and other uncovered microstructures the problem 
becomes more severe for the kind of covered microchannel structures used in the present 
invention, because the exciting and emitted radiation has to pass through plastic material. For 
30 transparent plastic material there is also a problem with "cross-talks" between the detection 
area/detection microcavities. Similar problems may also be at hand for spectroscopic methods in 
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which the radiation to be measured is created within the detection microcavity (for instance 
chemiluminescence, bioluminescence, etc.). 

A more recent problem relates to the fact that the present assignee recently has 
managed to control the liquid flow in microfluidic devices containing a plurality of microchannel 

5 structures in such a way that the inter-channel variation for a device with respect to flow 
becomes insignificant. This progress has enabled the assignee to quantify with a low inter-assay 
variation and a high sensitivity analytes, such as antigens, in the subfemtomole range in nl- 
volumes by carrying out the solid phase reaction of a heterogeneous sandwich immunoassay 
under flow conditions in small columns (nl-columns). This has raised the question about 

10 measuring the amount of an affinity complex such as an immune complex as a function of 
position along the flow direction in a column. See WO 02075312 (Gyros AB) and assignee's 
poster presented on September 17, 2001 at Proteomic Forum September 16-19, 2001, Munich, 
Germany. 

BRIEF SUMMARY OF THE INVENTION 

15 The present invention is directed to a method which meets the discussed problems. 

A first object is to reduce the electronic, chemical, optical and dust noise that is 
present in the stored measurement data of a search area, i.e. that is present in the stored image of 
a detection area. Electronic noise is introduced by the detector arrangement instrument. The 
optical noise arises from the design of the disc, the lid of the disc (if present), the detector head 

20 (e.g. lenses) etc. The chemical noise arise from chemicals and materials used in the manufacture 
of the disc and from undesired reactions including undesired adsorption related to the various 
reagents and solutions that have been used for placing the signal creating substance in the 
detection microcavity. Peak noise may derive from the manufacturing or the processing of the 
disc. Different kinds of dust, like particles, fibres, hair, etc may also introduce peak noise. The 

25 invention typically reduces noise-from steps related to chemical noise, e.g. disturbing 
fluorescence in the case the substance raw data image is based on fluorescence. 

A second object is to provide a method, software and arrangement for treating 
radiation data from one or more detection areas of the kind of microfluidic device described 
herein. 
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A third object is to provide arrangement, software and methods enabling accurate 
integration of radiation deriving from a desired substance as a function of subareas of individual 
detection areas. In particular this object aims at avoiding the problems discussed herein which 
the inventors have found may appear when measuring low amounts of substances with a high 
accuracy in microti uidic devices. 

These objects in particular apply to radiation data that have been collected while the 
detection area(s) is(are) moved in front of the detector head that is used for collection of the data, 
or vice versa, e.g. the detection area(s) is/are present in the surface of a circular disc that is spun 
in front of the detector head. 

The foregoing has outlined rather broadly the features and technical advantages of 
the present invention in order that the detailed description of the invention that follows may be 
better understood. Additional features and advantages of the invention will be described 
hereinafter which form the subject of the claims of the invention. It should be appreciated by 
those skilled in the art that the conception and specific embodiment disclosed may be readily 
utilized as a basis for modifying or designing other structures for carrying out the same purposes 
of the present invention. It should also be realized by those skilled in the art that such equivalent 
constructions do not depart from the spirit and scope of the invention as set forth in the appended 
claims. The novel features which are believed to be characteristic of the invention, both as to its 
organization and method of operation, together with further objects and advantages will be better 
understood from the following description when considered in connection with the 
accompanying figures. It is to be expressly understood, however, that each of the figures is 
provided for the purpose of illustration and description only and is not intended as a definition of 
the limits of the present invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present invention, reference is now made 
to the following descriptions taken in conjunction with the accompanying drawing, in which: 

Figure 1 illustrates a schematic view on an arrangement of the invention and its main 
parts. 

Figure 2 illustrates a detector head placed above a disc (microfluidic device) (cross- 
sectional view). 
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Figure 3 illustrates a set of microchannel structures that can be used in a circular 
disc. 

Figure 4 illustrates the search and detection areas. 

Figure 5 is a flow chart illustrating a first variant, a', of a method step a according to 
the invention. 

Figure 6 is a flowchart illustrating a second variant, a", of the method step a 
according to the invention. 

Figure 7 is a flowchart illustrating a second variant, a", of the method step a 
according to the invention. 

Figure 8 is a flowchart illustrating a method step % according to the invention. 
Figure 9 illustrates an edge curve. 

Figure 10 is a flowchart illustrating a method step 8 according to the invention. 

Figure 1 1 is a flowchart illustrating a method step 8 according to the invention. 

Figure 12 is a flowchart illustrating a method step § according to the invention. 

Figure 13 illustrates the search and detection areas and the surrounding areas A1-A3. 

DETAILED DESCRIPTION OF THE INVENTION 

It is readily apparent to one skilled in the art that various embodiments and 
modifications can be made to the invention disclosed in this Application without departing from 
the scope and spirit of the invention. 

As used herein, the use of the word "a" or "an" when used in conjunction with the 
term "comprising" in the sentences and/or the specification may mean "one," but it is also 
consistent with the meaning of "one or more," "at least one," and "one or more than one." 

As used herein, the term "rotating" refers to spinning. Yet further, the term 
"rotating" may also include, but is not limited to a step-wise rotation of a disc. 
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As used herein, the term "reagent" includes, but is not limited to an analyte. 

As used herein, the term "in a circular manner" refers to around the center of the disc 
(circumferential direction). 

The term "a plurality of microchannel structures" means two, three or more 
5 microchannel structures. Typically the term "plurality" means > 10, such > 50 or > 100 
microchannel structures. 

The terms "removing noise" and "reducing noise" means that the noise contribution 
in the measured data (values) is reduced to zero, or is negligible. The words "reduce" and 
"remove" will be used alternately in the text, but there is no difference in the meaning of the two 
1 0 words in this specification. They should be regarded as equivalent. 

As illustrated in figure 1, a detector arrangement is typically adapted for measuring 
radiation from a plurality of detection areas (31 7a.b.c. etc.. in figure 3) each of which is 
associated with a detection microcavity (306a,b,c, etc., in figure 3) in a microfluidic disc (101 in 
figure 1, 301 in figure 3). With reference to figure 1, the arrangement comprises: 

15 (a) a detector head (102) with a focal area, and a disc holder (105) which are linked to a 

means enabling for the detector head (102), i.e.. the focal area to transverse, the surface of the 
disc (101) in an essentially circular (means I) and/or radial manner (means II) when the disc is 
placed in the disc holder (105) 

(b) an angular aligning system (108.109) for recognizing the angular position of the part 
20 area which at a particular time is covered by the focal area, and 

(c) an optional radial aligning system (1 10,1 11) for recognizing the radial position of the 
part area which at a particular time is covered by the focal area, and 

(d) a controller (112), e.g, computer with software, which controls 

(i) means I and means II causing the focal area to transverse the detection areas 
25 (3 17a',b',c\ etc.) in an annular zone of the disc, and 
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(ii) the detector head (102) successively collects radiation in a preselected manner 
from individual subareas of essentially the same size as the focal area within at 
least one of the detection areas in said annular zone. 

In specific aspects, it is contemplated that the disc holder (105) is linked to at least 
one of the means for enabling the detector head and the disc to move relative to each other so 
that the focal area transverses the surface of the disc. Thus, the disc holder (105) is linked to 
either means I or means II or both means I and means II. 

Yet further, it is also understood by those of skill in the art that means I, means II, 
and the detector head can be linked directly or indirectly to the controller. 

Means I comprises three main variants with respect to how the disc and detector 

moves: 

(1) The disc is rotating around its axis of symmetry. In this variant, means I preferably 
comprises a spinner motor (103) and a shaft (104) carrying the disk holder (105). See figure 1 . 

(2) The detector head (focal area) is circularly moving around the axis of symmetry of the 
disc. In this variant, means I preferably comprises a spinner motor and a shaft with 
an arm that carries the detector head. This shaft typically has the same direction as 
the axis of symmetry of the disc when placed on the disc holder. 



(3) A combination of (1) and (2). 

Means II comprises three main variants with respect to how the disc and detector 
head moves: 

(1) The disc is laterally moving in front of the detector head. In this variant, means II 
preferably comprises a motor for laterally moving the disc holder. 

(2) The detector head (focal area) is laterally moving over the disc surface. In this 
variant, means II preferably comprises a linear frame (106) that carries the detector 
head, and a drive unit (107) for radial movement of the focal area, over a disc which 
is placed in the disc holder (105). See Figure 1. 

(3) A combination of (I) and (2). 
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The arrangement is illustrated in figure 1 . 

I. Means I and the Angular Aligning System for the Detector Head 

In a typical variant, the detector head (102) and the motor (103) (e.g., a spinner) with 
a rotatable shaft (104) carrying a disc holder (105) are supported on a frame structure (113). The 
5 motor (103) controls the rotating speed that can be varied, e.g., within an interval between 0 - 
15,000 rpm, such as within an interval between 60-5,000 rpm. The rotation of the disc may be 
stepwise. When the disc is rotating, the focal area of the detector head will successively scan 
angularly adjacent part areas of the disc surface. Means I in figure 1 is according to variant (1) 
above and comprises the motor (103) and the shaft (104). 

! 0 The disc holder (1 05) is preferably a plate on which the disc can be placed. The disc 

holder could also be a device that holds the disc at its periphery. In order to reduce wobbling of 
the disc (if the disc is skewed), the side (114) of the plate facing the disc may comprise a system 
of evenly distributed uncovered shallow grooves or channel openings that are connected to a 
vacuum system by which the disc can be sucked to the plate. See for instance PCT/SE02/01682 

15 (Gyros AB). 

Detection areas may be present in one or both sides of the disc. If the disc holder is in 
the form of a plate as illustrated in figure 1 and the microfluidic disc has detection areas on the 
side opposing, it is important to secure the plate so as not to disturb the collection of radiation. 
The plate may thus have a smaller diameter than the disc with the detector areas being located in 
20 an annular zone that is not covered by the plate. Alternatively, the plate may be in a material that 
is translucent for the radiation utilized for the measurement, at least at the detector areas. 

The angular aligning system may comprise: 

(1) a device that will enable the determination of when a predetermined angular position 
of a disc placed on the disc holder is in front of the objective of the detector head 

25 (102) (i.e., covered by the focal area), and/or 

(2) a home position mark detector (108) which is able to detect when a home position 
mark (305) on a rotating disc (101) placed on the disc holder (105) is passing. 

A home position mark (305 in figure 3) is preferably placed in an outer 
circumferential zone outside the detection areas or in some other position, which can be detected 
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with high accuracy. The position coordinates of each specific spot of the surface of the disc is 
given as the angular position relative to the home position mark and as the radial position 
relative to the circumference or axis of symmetry or relative to any other arbitrary fixed position 
on the disc. 

A home position mark detector (108) typically has a fixed position outside the disc, 
for instance on the frame structure (113). 

An accurate and preferred alternative for determining when a predetermined angular 
position is in front of the objective is to include an encoder that progressively gives the angular 
distance from the home position mark while the disc is rotating. This kind of encoder (109) is 
typically associated with means I, for instance the motor (103), the shaft (104) or the disc holder 
(105). By associating the encoder directly with the disc (101) it is likely that the most accurate 
determination will be accomplished. The encoder typically divides each revolution of the shaft 
into a large number of grades, for instance > 5 000, such as > 10 000 or > 20 000 or > 30 000. A 
simple but less accurate alternative is to include calculating means that calculates the time 
1 5 needed from a preset rotation speed and the angular distance between the predetermined position 
and the home position mark (i.e., from the preset rotation speed and the angular position co- 
ordinate). This kind of calculating means may be associated with the controller. 

The angular aligning system should be able to give the angular position coordinate 
for the part of the disc which is in front of the objective with an accuracy of ± 1°, such as within 
20 ± 0.1° or within ± 0.01° (provided there are 360° per revolution). The exact accuracy needed will 
depend on the size of the disc, radial position of the detection area, the required sensitivity, size 
of detection area, etc. 

TT Means II and the Radial Aligning System for the Detector Head 

H. Mea ^ e U detector head (1Q2) * guided on a linear frame (106) that may be the upper part 

95 of the frame structure (113) for linear displacement and positioning in a first plane PI, 
transversely through the centrales CL of shaft (104) and running in a radial direction thereto. 
The linear frame (106) prohibits uncontrolled movement of the detector head in any other 
direction relative to this linear displacement. The drive unit (107) for this displacement may be 
in the form of a translation^ responder for incrementally changing the position of the detector 

30 head (102) in the first plane PI (radial movement) and for enabling scanning of radially adjacent 
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subareas of a microfluidic disc device placed in the disc holder (105). Means II comprises the 
linear frame (106) and the drive unit (107) and is in figure I according to variant (2). 

The drive unit (107) is associated with a unit (110) for determining the linear 
displacement and thus also the radial position coordinate of the part of the disc which is in front 

5 of the objective of the detector head (102). This unit (110) may be in the form of an encoder that 
gives a translational position and movement of the focal area (objective of the detector head) in 
relation to a translational home position (1 1 1) on the responder. This home position in turn may 
be associated with a unique radial position in the disc. The measuring unit (110) should be able 
to translate a translational position and movement of the detector head (focal area) into a radial 

10 position coordinate of the disc used with a high accuracy, typically within ± 10 um such as ± 1 
urn or ± 0.1 (-im. 

The drive unit (107) and the vertical height of plane PI may be adjustable for 
focusing purposes. 

15 11 L C ° n Control means, for instance electronic and programmable control means 
(schematically illustrated by reference numeral (1 12)) with operator's interface and software, not 
further disclosed, may be assigned to the detector arrangement among others for 

a) recognizing one or more pairs of start/stop-positions (angular and/or radial) for 
irradiating if the detection principle utilized requires irradiation and/or for collecting 

20 radiation, 

b) identifying individual subareas in detection areas or elsewhere in the surface of the 

disc, 

c) controlling the simultaneous rotating of the disc and the incremental radial 
displacement of the detector head (102), 

25 d) collecting radiation data from the detection areas/detection microcavities. 

e) treatment and presentation of the collected data, and/or 

f) determining the time at which a particular angular position is in front of the objective 

of the detector head from the rotational speed. 
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g) storing and reading collected data in a data storage 130, 

h) storing a computer program/computer program product, i.e. a computer software 
application in a computer software application and computer program storage 120, 
and reading and accessing software applications from the same. 

5 The present invention may be implemented as a computer program product, i.e. a 

computer software application, directly loadable into an internal memory storage of a processing 
unit within the controller 120. Said computer program product comprises the software code 
means for performing the steps of the invented methods. 

10 Further, the present invention relates to a computer program product stored on a 

computer usable medium, comprising readable program for causing a processing unit in a 
computer means, such as the controller 1 12, to control an execution of the steps of the invented 
methods. 

15 The computer usable medium is a record medium, a computer memory, a Read-Only 

Memory or an electrical carrier signal. 

Different parts of the arrangement may communicate (115) with the controller (1 12). 
The controller will in the preferred variants instruct the detector head to successively collect 

20 radiation from distinct and preselected parts of the surface of the disc. Typically the controller is 
programmed to start collecting radiation at a position, primarily an angular and/or a radial 
position, which is prior to an intended detection area, and to stop collecting at a position, which 
is after the same detection area. The start and stop positions define a search area that at least 
should encompass the detection area or the relevant part of the detection area in which there is 

25 radiation from the substance. See further figure 4. If the radiation requires that the substance is 
irradiated, which is the case if-fluorescence is measured, the control means also defines the 
settings for the start and stop positions for irradiation. These latter settings are typically 
essentially the same as for collecting radiation. 

The start and stop signals for collecting radiation is preferably directly linked to the 
30 angular positions in the disc at which collection is to start and end, respectively. This also 
includes that due account is taken for delays that may be inherent in the system or preset, i.e., the 
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start and stop signals may have to be initiated before the focal area is positioned in front of the 
start and stop position, respectively. If the angular aligning system comprises an encoder, the 
encoder signals corresponding to a start position and a stop position are used to define the time 
' -period during which radiation is. to be collected. In an alternative, the start and stop for collecting 
5 radiation is linked to a preset rotating speed, i.e., the controller calculates from a preset rotating 
speed the time at which the start and stop position should be in front of the objective. 

The controller may be programmed to change the radial position of the detector head 
(focal area) after a predetermined number of revolutions of the disc, for instance after 1. 2 or 
more revolutions with preference for 1. 

10 The controller may be capable of changing the radial position of the objective (focal 

area) during a revolution of the disc. For these variants, one can envisage that radiation is 
collected for all relevant subareas at a common angular position before the disc is rotated (in a 
single step) to a subsequent angular position. In an alternative variant the objective (focal area) 
is transversing the disc surface in a spiral-like manner, i.e., the radial position is changed 

1 5 successively during a revolution. 

In the preferred variants, the collected radiation data is stored in a data storage 130 in 
a form that is retrievable for each individual subarea, for instance in the control unit. This means 
that after collecting of radiation, it will be possible to represent the collected data as a 3-D image 
of each search area and/or detection area showing the amount of radiation from each individual 
20 subarea. In the case of overlapping subareas, the proper treatment of the data takes into account 
the overlapping effect and creates a true image of the radiation associated with different parts 
(subareas) of a detection area. 

IV. Detection Principles and the Detector Head 

The microfluidic devices used in the present invention typically require detection of 
25 very low absolute amounts or concentrations of substances in the detection microcavity. It is 
therefore imperative in many variants of the invention that the detection principle shall enable 
detection and quantification of substance amounts that are < 10" 12 mole per detection 
area/detection microcavity, such as < 10"' 5 mole or < 10' 18 mole per detection area/detection 
microcavity. 
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Typical detection principles may include, for example, spectrometric detection, such 
as radioactivity, absorption of radiation (e.g. light), fluorescence, chemiluminescence, 
bioluminescence, scattered light etc. For variants in which radiation requires irradiation a change 
in wave-length(s), polarization, life-time, scattering, intensity, etc., between irradiation and 
radiation may be measured as a function of the presence of the substance of interest in the 
detection microcavity. Typical examples include, but are not limited to luminescence and 
fluorescence principles, with laser induced fluorescence (LIF) being preferred. 

Details about the detector head are given in copending applications PCT/SE02/01678 
and corresponding US application S.N. 10/062,258. 

The size of the focal area of the detector head typically is less than the size of a 
search and/or a detection area. The width of a focal area in one direction (direction 1) is typically 
essentially < 1/5, such as of the < 1/10, of the corresponding width of a search and/or detection 
area and in a perpendicular direction (direction 2) within the same limits or larger. A suitable 
detector head in the form of a pick-up head (200) is illustrated in figure 2 and comprises a laser 
source (201) whose beam is reflected on a dichroic mirror (202) and focused through an 
objective (205) to a part of a detection microcavity (203) positioned in front of the head. The epi- 
fluorescent light is passed through the dichroic mirror and through a band-pass filter (206), 
selective for the flourochrome at hand and is finally focused onto the entrance of a photo 
multiplicator tube (PMT) (207) by means of an aspheric lens (208). Pin-holes (209 and 210) are 
positioned between the entrance of the laser beam and the dichroic mirror, and in front of the 
PMT (207). The size and position of the pin-hole (210) are adapted so that the focal area of the 
laser beam is inside the detection microcavity. The size of pin-hole (209) is adapted so that 
preferentially emission light emanating from the focal area is passed into the PMT (207). A 
mirror spot on a glass disc or the like may replace the dichroic mirror. 



V. Microfluidic Device 

The microfluidic device used in the various aspects of the invention comprises a 
plurality of microchannel structures in which aliquots of liquids are transported and/or processed. 
The devices typically are disc-shaped. The structures are covered in the sense that their interior is 
in contact with ambient atmosphere primarily via separate inlet and/or outlet openings and/or 
vents. Each microchannel structure comprises one or more detection microcavities and possibly 
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also one or more reaction microcavities, and microconduits connecting these parts with each 
other. A reaction microcavity may coincide with a detection microcavity. The result of the 
processing in a microchannel structure is measured as radiation from a detection area which is 
directly or indirectly associated with a detection microcavity. 

5 A disc typically has an axis of symmetry (C„) perpendicular to the disc plane where n 

is an integer > 2, 3, 4, 5 5 preferably co (C„). The disc may thus be rectangular or have a circular 
shape. 

Different principles may be utilized for transporting the liquid aliquots within the 
microfluidic device/microchannel structures. Thus inertia force may be used, for instance by 
10 spinning the disc (centrifugal force). Other forces are capillary forces, electrokinetic forces, 
hydrodynamic forces etc. 

Microfluidic devices that have an axis of symmetry and are intended for rotation may 
have a home position mark as discussed above. 

The terms "microchannel", "microconduit", etc., contemplate that a channel structure 
15 comprises one or more cavities and/or channels/conduits that have a cross-sectional dimension 
that is < 10 3 urn, preferably < 10 2 urn. The lower limit is typically significantly larger than the 
size of the largest reagents and constituents of aliquots that are to pass through a microchannel. 
The volumes of microcavities/microchambers are typically < 1000 nl. such as < 500 nl or < 100 
nl or < 50 nl or < 25 nl, which in particular applies to the detection microcavities. 
20 Chambers/cavities directly connected to inlet ports for liquids may be considerably larger, e.g.. 
microchambers/microcavities intended for application of sample and/or washing liquids. 
Microformat means that one, two, three or more liquid aliquots that are transported within the 
device have a volume in the ul-range, i.e., < 1000 ul such as < 100 ul or < 50 ul including but 
not limited to the nl-range (nanoformat). such as < 1000 nl or < 500 nl or < 100 nl or < 50 nl. 

25 The microfluidic device may be made from different materials, such as plastics, 

glass, silicone polymers, etc. The detection area should be transparent/translucent for the 
detection principle utilized by the detector. From the manufacturing point of view plastic 
material is many times preferred because the costs for the material are normally low and mass 
production can easily be done, for instance by replication. Typical manufacturing processes 

30 involving plastic material are replication by embossing, moulding, etc.. followed by attaching a 
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top lid covering the open microchannel structures so obtained. See for instance WO 9116966 
(Pharmacia Biotech AB, Ohman & Ekstrom). However, plastic materials may interfere with 
several sensitive detection principles. Their high auto-fluorescence is disadvantageous for 
normal fluorescence techniques in case low absolute amounts of fluorescent substances are to be 
5 measured. This points to the fact that it is important to match the material in the device with the 
detection principle used. At the priority date of this invention the preferred disc material is 
plastic material, such as polycarbonates and plastic material based on monomers which consist 
of a polymerisable carbon-carbon double or triple bonds and saturated branched straight or cyclic 
alkyl and/or alkylene groups. Typical examples are Zeonex™ and Zeonor IM from Nippon Zeon, 
10 Japan, with preference for the latter. See for instance WO 0056808 (Gyros AB, Larsson, Ocklind 
and Derand) which is hereby incorporated by reference. In this context silicone polymers such as 
poly dimethyl siloxane (PDMS) and the like are not considered to be plastic material. 

It is known that black plastic material, for instance containing graphite powder or 
carbon black, absorbs light and therefore has a low auto-fluorescence. Transport of light within 
15 black plastic material is prevented. Black plastic material will be very efficient for micro fluidic 
devices when fluorescence and luminiscence measurements are relied upon. Black plastic 
material should in this case be avoided in the detection areas. 

In case a lid is needed as for most discs obtained by replication, plastic lids of 
different origin may be used, for instance Melinex™ 12 PET and Melinex™ 17 OPP (Du Pont, 
20 U.S.A.), etc. 

From the auto-fluorescent point of view an optimal combination of transparent 
plastic material appears to be Zeonor !M for replication and Melinex IM 17 OPP as the lid. This 
combination of material is likely to be useful for excitation wavelengths in the interval 480-650 
nm, 

25 The plurality of dejection microcavities and the corresponding search and/or 

detection areas are preferably arranged in subgroups such that all members of a subgroup are 
positioned at the same radial distance and/or at the same angular position and/or have equal 
length and/or cross-sectional dimensions. Within each subgroup there may be at least two, three 
or more detection microcavities (detection areas), such as > 10 or > 25 or > 50 detection 

30 microcavities (detection areas). 
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A detection area in the inventive arrangement typically has a size within the range of 
lxlO 2 - 2xl0 6 urn 2 , such as lxlO 3 - 10 5 um 2 . Their length and/or breadth are typically within the 
range of 0.5x10 - 5xl0 4 um, such as 1x10 -10 4 urn. 

Design microchannel structures (301a,b.c. etc.) that are adapted to a circular disc are 
illustrated in figure 3. The structures are linked together by a common distribution channel (302) 
and a common waste channel (303). The orientation of the microchannel structures around a 
common axis of symmetry is apparent. The circumference (304) of the disc has a home position 
mark (305). Each of the combined reaction/detection microcavities (306a,b,c, etc.) is 
communicating in the downstream direction with the common waste channel (303) and in the 
upstream direction via separate connections with the common distribution channel (302) and 
separate volume measuring units (307a,b.c, etc.). A surface detection area (317a,b,c, etc.) is 
associated with each detection microcavity. The common distribution channel (302) carries at 
one of its ends and at an intermediary position inlet ports (308 and 309. respectively). Another 
kind of inlet port (310) is located at each volume measuring unit (307a,b,c, etc.). Each 
microchannel structure (301a.b.c. etc.) also has an outlet to the common waste channel (303) and 
an outlet port (318) at the remaining end of the common distribution channel (302). An inlet vent 
(311) to ambient atmosphere is connected to the common distribution channel via a common 
venting channel (31-2). Other vents (313 and 314) to ambient atmosphere are placed in the 
common waste channel (303) and in the connecting microconduit between each 
reaction/detection microcavity (306a,b,c. etc.) and the common waste channel (303). Appropriate 
valvings are positioned at 315 and 316 in each microchannel structure (301a,b,c, etc.). 

The diameter of the preferred circular discs is about the same as conventional CDs 
but may be larger or smaller, for instance up to 300 % or larger and down to 10 % or less. 

VI. Processes to be Preformed within the Microfluidic Device 

The processes that are carried out within the individual microchannel structures 
comprise assay protocols, organo-^hemical or biochemical synthesis protocols, etc. Typically the 
protocols comprise introduction of one or more liquid aliquots containing the necessary 
reagents/reactants into a microchannel structure. In the case of assay protocols, one of the 
aliquots is a sample which is uncharacterized with respect to at least one feature, e.g.. type, form 
and/or amount of an analyte. 
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The processes comprises that the substance, which is associated with the radiation to 
be collected, is formed and/or retained in the detection microcavity under static conditions or 
under flow conditions. The reaction system for retaining may be homogeneous or heterogeneous, 
i.e., with or without the desired substance being partitioned between a liquid phase and a solid 
phase. In case of flow conditions and microfluidic devices in the form of discs, the flow direction 
in the detection microcavity may be towards the circumference (outwards) or towards the center 
of the disc (inwards), or essentially parallel with the circumference of the disc. Also other 
directions may be utilized. 

Process protocols may utilize specific reactions between reactants having mutual 
affinity to each other leading to a (a) formation of an affinity complex that is immobilized to a 
solid phase in a detection microcavity or (b) one or more other reaction products that may be 
soluble or insoluble in the detection microcavity. By properly selecting the reaction conditions 
including selection of reactants, it can often be arranged so that the product obtained and/or a 
reagent in excess are detectable with a signal that can be (a) measured from the above-mentioned 
detection areas and (b) related to one or more features of a starting liquid aliquot introduced into 
a microchannel structure. Typical such features are kind, form and amount including activity. 
etc., of a particular reactant (e.g. an analyte) including for instance an affinity reactant such as an 
enzyme etc. The term "can be related to one or more features" includes also the determination of 
the manner in which reaction variables such as pH. ionic strength, detergents, etc., might 
influence the reaction used for forming the reaction product. Typically, one makes use of 
detection principles based on radioactivity, fluorescence, chemiluminescence, bioluminescence, 
enzymatic activity, chromogens, light scattering (turbidometry), etc., for instance by utilizing a 
reactant that carries a group providing detectability, either by being detectable as such or by 
being transformable to a detectable group. Typically the utilized process protocol means that a 
detectable reactant is incorporated into a complex or into some other reaction product. See for 
instance WO 02075312 (Gyros AB). Detectable products, reagents, etc., that can be formed 
and/or retained and measured iq a microcavity via a detection area are collectively called 
"substance" in other parts of this specification. 

Typical reactants in this context include members of affinity pairs such as (a) 
antigen/hapten and the corresponding antibody including its antibody active fragments, (b) lectin 
and the corresponding carbohydrate structure, (c) native ligands and the corresponding receptors, 
(d) complementary nucleic acids including synthetic variants such as synthetic oligonucleotides. 
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(e) Ig(Fc)-binding proteins and Protein A, Protein G and other Ig(Fc)-receptors, (e) ion pairs of 
opposite charges, enzyme and the substrate, inhibitor, cofacton coenzyme etc., that can bind to 
the enzyme, etc. Synthetic variants more or less mimicking a native affinity interaction are also 
included. 

5 One aspect of the invention is a method and corresponding arrangement, for 

determining the amount of a substance in a detection microcavity of a microfluidic device and 
comprises collecting radiation associated with the substance from a detection area associated 
with the detection microcavity. The method is characterized in comprising the steps of: 

a) providing 

10 (i) a microfluidic device, e.g., in the form of a disc, comprising 

A) a plurality of microchannel structures, each of which has an inlet port, 
a detection microcavity and a microconduit connecting the inlet port 
with the detection microcavity, and 

B) a plurality of detection areas and associated search areas as discussed 
15 elsewhere in this specification, each of the detection areas being (1) 

associated with one of said detection microcavities, (2) present in the 
surface of said device and (3) translucent/transparent for said 
radiation, and 

ii) a detector arrangement which is capable of collecting radiation from 
20 individual subareas of each of said search and detection areas; 

b) processing one or more liquid aliquots in at least one of said plurality of microchannel 
structures so that the substance is retained in the detection microcavity of each of said 
at least one of said plurality of microchannel structures, provided that said substance 
and/or one or more reagents that are necessary for the substance to be retained in a 

25 detection microcavity are present in at least one of said one or more aliquots; 

c) scanning the search and/or detection areas associated with the detection microcavities 
that are part of microchannel structures in which step (b) has been carried out to obtain 
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radiation from individual subareas of each scanned search or detection area, said 
scanning being performed by the use of said detector arrangement; 

d) integrating radiation as a function of the subareas of each scanned search area or 
detection area to obtain the amount of radiation from each search area or detection 

5 area; 

e) characterizing for each of the amounts obtained in step (d) a reaction variable that has 

been included in the process protocol used for each microchannel structure. 

Further details about characterization of reaction variables are discussed above and in 
WO 02075312 (Gyros AB). The characterization includes e.g., that the amount of the substance 

10 in absolute or relative terms in each of the detection microcavities is determined from each of the 
amounts obtained in step d). The substance in this context may be the substance from which the 
radiation is directly derived or a substance, such as an analyte, the amount of which in one of the 
liquid aliquots used in step (b) is related to the amount of radiation determined in step (d). Since 
it is the result of the reaction that is measured the analyte may or may not be present in a 

1 5 detection microcavity after step (b). 

Step (c) comprises collection of radiation intensities which refer to changes in the 
radiation due to the amount of a particular substance in the detection microcavity. This includes 
true intensity values but also other changes, such as changes in refraction, in scattering, in the 
plane of polarization and life time of fluorescent light (including half life, e.g.), etc. 

20 The microfluidic disc and steps b) and c) are illustrated elsewhere in this 

specification. In addition to circular scanning, e.g., by rotating a disc, step (c) also comprises 
non-circular scanning for instance scanning by lateral movement of a linear detector head 
comprising one or more rows of detector elements over the detection areas. Also imaging by a 
CCD camera is included. The microfluidic disc may be circular or have some other geometric 

25 form, for instance triangular, rectangular, etc., including also irregular forms. 

Step (d) means integration over each detection area, Le., primarily over subareas 
which have radiation values that deviate from the values obtained for surface parts of the device 
that surrounds the detection area. Alternatively one may exclude all or selected parts of a 
detection area which corresponds to parts of a detector microcavity in which the presence of the 
30 substance is insignificant. The integrating step includes the substeps of (a) finding a start and/or 
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stop position at edges of the detection area, e.g., at an end corresponding to the inlet end of the 
detection microcavity, and (b) the factual integrating. In one variant substep (a) is carried out by 
determining the inflection point for the amount of radiation per subarea versus position along the 
detection area. Due the fact that the delineating part of the detection area may be curved, at the 
upstream end, the invention also suggests that the integration should account for curvatures in 
the circumference of a detection area. In a preferred variant, substep (a) comprises determining a 
threshold that segments detection area pixels from the background, e.g., if the detection 
microcavity is filled with a particle bed a threshold that segments the particle bed pixels from 
background pixels. This can be done with optimal thresholding or determining median or mean 
background or any other way of determining background (for "optimal thresholding" see for 
instance reference [1]: "Digital Image Processing", 2nd edition. Editors: Gonzales R C et aL 
page 354). In this case substep (b) will mean integrating a selection of pixels, Le. t those pixels 
which have radiation values above the threshold and belong to the main group (detection area) 
and excluding noise pixels above threshold that do not belong to the main group of connected 
pixels. The integrating typically starts from pixels corresponding to one end of the detection 
area, e.g., the inlet end. 

Between the preferred variant of substeps (a) and (b), there are preferably additional 
substeps for refining the method, such as 

i) Creating a binary image from the calculated threshold. 

ii) Labelling the high binary pixels into different groups (= labeling the image). Each 
group will consist of pixels that are bordering to each other (close to each other). The 
binary high pixels from the detection area will define the main group. 

iii) Optionally removing all binary high pixels that do not belong to the main group. 

The factual integrating (main step d) will in this case mean integrating the radiation 
values for the binary high pixels of the main group. 

The inventors have found that when working with very small volumes and amounts, 
the material from which a device is manufactured as well as the pretreatment procedures may 
introduce radiation artifacts in form of peak noise that is comparable to the radiation coming 
from subareas outside the peak. Therefore step (d) may also comprise substeps for removing 
peak noise. This typically means a first substep in which the deviating radiation (peaks) are 
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made more apparent. One way of doing this is la Place filtering, point detection and the like. See 
for instance see reference [1]: "Digital Image Processing", 2nd edition, Editors: Gonzales R C et 
al, pages 333, 334 and 339 (1987). In a subsequent substep the width and the position of each 
peak is calculated for instance by including edge detection or edge linking based on a local area 
5 etc.. See for instance reference [1]: "Digital Image Processing". 2nd edition. Editors: Gonzales R 
C et al, pages 334 and 344 (1987). In the next substep the peak noise is removed by interpolating 
from surrounding subareas. An alternative way to the whole process of removing peaks is 
morphological opening, un-linear filtering operating in local histogram domains etc. These 
substeps for removing peak noise is performed prior to the factual integration. 

10 Step (e) is conventional and typically includes that the integrated value is compared 

with the value for one or more standards. A standard value is typically the integrated value for a 
known amount of a standard substance, which in most cases is the same as the substance under 
investigation. 

The scanning step (c) and integrating step (d) may in certain innovative variants of 
15 one aspect of the invention be optional. Thus, these steps are preferably included when the 
desired substance is unevenly distributed within the detection microcavity. In a typically case, 
this may happen when the substance is retained within the detection microcavity under flow 
conditions, for instance by being captured to a solid phase introduced into the microcavity prior 
to the fluorescent substance. 

20 In the that case the substance is homogenously distributed within the detection 

microcavity, the scanning step (c) and integrating step (d) may be replaced with collecting the 
radiation intensities for selected subareas of a detection area and letting these intensities 
represent the total amount of radiation from the detection area, for instance as a mean or 
maximum value. Step (e) can then be carried out on these values in the same manner as for 

25 values obtained by scanning and integrating. This way of performing the method is also a part of 
the present invention. * 

Homogenous distribution of the substance in the detection microcavity typically is at 
hand in case the substance is present in the detection microcavity in an equilibrated solution 
and/or when a reaction is going in the microcavity between homogeneously distributed reactants 
30 either forming or producing the substance that is to be detected. 
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The innovative method of the aspect discussed in the preceding paragraphs includes 
that steps (d) and/or (e) are performed in close connection to the preceding steps [(a)-(c) or (a)~ 
(d), respectively] and/or that the data from the scanning and/or the integrating have been 
obtained at an earlier time, for instance at a different geographical location, and/or by different 
individuals. 

In the described method for determining an amount of a substance in a detection 
microcavity of a microfluidic device, it may be necessary to reduce the noise contribution in the 
detected radiation intensity from scanning step (c) in the integration step (d). 

A preferred embodiment of the invention is now described. 

Method for reducing noise in a substance raw data image containing noise. 

An algorithm, or rather an automated method, for that is detectable by use of a particle bed is 
presented hereinafter. This algorithm may be started for each scanned search area/detection 
area/particle bed of interest. 

The invention will now be further detailed with respect to a particular example comprising a 
microfluidic device (see figure 3) in a disc format, comparable to Compact Disc (CD), having an 
axis of symmetry perpendicular to the disc plane and with a radial flow direction from the axis of 
symmetry for each of the detection microcavities (see figure 3: reference number 306 a,b,c . . . ) 
that comprises a bed of particles. This kind of microfluidic device will henceforth be called CD. 
The signal detected and measured may be fluorescence from an antibody fluorochrom conjugate 
(antibody-fluorochrom). Each bed defines a detection area (see figure 3: reference number 317 
a,b,c ... ) in the surface of the disc. See for instance WO 02075312 (Gyros AB). The detection 
area is part of a search area comprising also a defined minor zone around the detection area. The 
principles outlined in this example may be adapted also to other devices. The principles are also 
applicable to other kinds of radiation as discussed elsewhere in this text. 

Generation of a substance raw data image 

As described in section III, Controller, of this description, a controller means (fig. 
1: ref. No. 1 12) is programmed to control the scanning and detection. 
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Each search area is scanned by a detector head that is capable of detecting a 
detectable physical feature deriving from a substance present in the microcavity, like radiation 
caused by radioactivity, fluorescence, chemiluminescene, bioluminescence, enzymatic activity, 
scattering of light, absorption of light, reflection of light etc. The measured data are stored in a 
5 data storage, an image data storage, connected to or connectable to an processing unit, which 
may either be comprised in the controller means, or as in another possible embodiment, 
incorporated in a separate computer system. Each search area is identified by an identifier, which 
may be used as address. The identifier makes it possible for the system and controller means to 
locate each single search area and its detection area on a CD and link the corresponding stored 
10 measured data in the image storage. 

Now with reference to figure 4. each search area 402 has a pre-defined radial start 
position, Xstart, a pre-defined angular start position, YstarL a pre-defined angular stop position. 
Ystop, and a pre-defined radial stop position, Xstop. The start of the collecting of data is started 
at a position, primarily an angular and/or a radial position that is prior to an intended detection 

15 area, and to end the collecting at a position that is after the same detection area. The search area 
will therefore cover an outer area 406 surrounding the detection area 404. The search area 402 is 
scanned in radial direction by the detector disc lap by disc lap. The radial start position Xstart is 
the number of disc laps/revolutions, counted from a spinning start position, Spinstart (not 
shown). When the scanning of a search area starts, and the radial stop position is the number of 

20 laps, counted from the spinning start position, when the last angular scanning of the search area 
is to be performed. The scanning of the search area in the angular direction starts during the lap, 
which is defined by the radial start position, at the angular start position in the angular direction, 
which is defined by the spinning direction of the disc and the detector position in relation to a 
Home Mark on the CD. The scanning in angular direction is stopped at the angular stop position. 

25 Figure 4 is schematic in the sense that the search area of a rotatable device typically 

is a) slightly narrowing when going from an outer position inwards towards the centre of the CD, 
and b) slightly arc-shaped. Compare for instance that the Ystart position for each lap are always 
at the same angular position (which also applies to the Ystop position). 

During a lap, from the angular start position to the angular stop position, a stream 
30 of collected data from that search area will generate a corresponding stream of measured data 

that are stored into the storage. This is repeated lap by lap from the radial start value to the radial 
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stop value. The processing unit of the controller 120 links measured data to the correct identified 
detection area, and there will be no mixing of measured data from different search and/or 
detection areas. Each scanned search area has a corresponding set of collected measured data. 
The processing unit will organize the data in such an order, position by position in each set of 
5 measured data, and it is therefore possible for the image processor to locate a certain measured 
value in the set. The set of measured data has therefore similarities to an electronic and/or 
digitally stored image. An image is defined in very small elements, called pixels. Each pixel 
corresponds to a position of a measurement area in the search area. 

When the scanning, detecting and storing are finished, the measured data of a 

10 search area is possible to plot by use of printer/plotter software program used by the processing 
unit in a three dimensional diagram, wherein the horizontal plane is defined by the two 
dimensions x and y, in this case the radial direction, x-axis, and the angular direction, y-axis, 
which is perpendicular to the x-axis. Positive x-direction is defined as the main flow direction 
(downstream) through the detection area, from the inlet 408 of the detection area 404 to the 

15 outlet of the same. Positive y-direction is defined as the scanning direction which is 

orthogonal/perpendicular to the radius (or flow direction). The third dimension, z-axis. 
perpendicular to the horizontal plane defined by the x- and y-axis, is the magnitude of the 
measured value, e.g. measured value of the fluorescence intensity. The measured data defines a 
surface of measured data, having peaks and valleys. As mentioned above, the three dimensional 

20 figure has essential similarities with an image having an intensity value defined for each pixel, 
and the surface of measured data is therefore called an image of a search area/detection 
area/particle bed comprising the substance from which the measured signal derives. In the 
description of this application, an image is the stored detector measurement data belonging to a 
search area/detection area/particle bed. A pixel is here defined as a discrete measurement area 

25 within the search area/detection area and each pixel may have a measured intensity value, a pixel 
value. The resolution of the image is depending of the number of discrete measurement points in 
the y-direction and the number of teps in the x-direction within the search area. The resolution 
increases with increasing density of measure points and scanning laps, and consequently, the 
area of a pixel will decrease. If the electromagnetic-emitting source is a laser, the possible 

30 minimum pixel area is limited by the minimum cross-section of the laser beam with the search 
area. 
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During the scanning of the search areas the pixel intensities are determined by the 
detection process and stored as digital pixel data representing raw data images of the search 
areas. 

The space between the surface of measured data for the search area/detection area 
and the base plane, defined by the two axes in the base plane, constitutes a volume. The whole o 
a part of the volume may represent a wanted measure for a chemical process, i.e. a reaction 
variable. This value is possible to calculate by using a mathematical method, e. g. addition or 
integration of measurement data over the whole area of the detection area, or one sub-area 410 
(fig. 4) or a number of sub-areas of the detection area. Sometimes the object is to calculate an 
amount for an subarea, here called a normalized area, i.e. a predefined amount of pixels, inside 
the detection area. This is done to minimize the variance between detected structures. The 
variance depends among other things on background radiation, such as fluorescence, from the 
detection area. 

In an ideal image of a detection surface, the area is smooth, because of the absence 
of noisy background contribution in the measured raw data signals from the detector. This image 
includes only measured data from a substance signals, i.e. the substance contribution to the 
detected and measured raw data signal. In reality, a substance raw data image comprises 
measured substance data and noise data. The noise cause irregularities, like noise peaks, in the 
surface, and evidently, a calculation of the volume in an area having one or many noise peaks 
will not result in a true measure for a chemical process. 

It is therefore necessary to reduce the noise contribution in a substance raw data 
image as much as possible before calculation of the measure. 

The present invented method provides a method for reducing noise contribution in 
a substance raw data image derived from a detection area of a microfluidic device of the kind 
described above. The method comprises a number of main-steps for reducing different kind of 
disturbances and noise within an image. 

The method removes noise from substance raw data information in substance raw 
data images. Each image is associated to one detection microcavity and one substance that is 
present in the detection microcavity of a microfluidic structure of a microfluidic device. The 
substance raw data information is obtained by use of a detector arrangement capable of scanning 
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a search area associated to a detection area, which is associated to the detection microcavity. The 
detector arrangement is capable of detecting at least one of a number detectable physical features 
discussed above, e.g. fluorescence, luminescence, etc. emitting radiation from the search area, 
and generate output that is proportional to the magnitude of detectable activity. The output is 
converted to measure data storable as substance raw data in a set in a data storage. The method 
may comprise following steps in any order for removing different noise signals contributions 
step-by-step in any order in a substance raw data image: 

- reducing background radiation (step a); 

- reducing peak disturbance (step (3); 

- locating the detection area within the search area (step x); 

- moving binary artefacts (step 8); 

- removing unwanted areas of the detection area (step s); 

- applying default detection area in noisy images (step (j>). 

Each of the steps a - (|) comprises a number of sub-steps. Steps a - <j) may be regarded as 
an algorithm, but is not necessary that all of the steps a - <j> are carried out for all images. The 
order and kind of steps is typically preselected and depends among others of the sensitivity and 
accuracy required in each individual step. In general, inclusion of more steps is required when 
high sensitivity is desired. 

The different steps will now be described in more detail. 



Step a: - Reducing background radiation and construction of background 
processed images. 

The object is to remove the electronic, optical and dust noise that may be visible i 
the image. This part of the invented method mainly comprises two variants. 

In the first variant, a background radiation image is used. In figure 5, the first 
variant is presented schematically in a flowchart. Following steps may be used: 

- if available, accessing a background radiation image or. if not available, 
generating a background radiation image (step a' 1); 
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- if available, accessing corresponding substance raw data image or, if not 
available, generating substance raw data image (step a'2); 

correlating the two images and generating a correlation image (step a' 3); 

- locate the positions of the largest pixel values (step cc'4); 

- positioning the images correctly, subtracting corresponding pixel values and 
storing the generated background processed image(step a' 5). 

This first variant is now discussed in more detail. 

In the first variant the First step, step a' 1, is to scan the disc and measure the 
radiation, e.g. the intensity of the fluorescence, without presence of the substance giving raise to 
the measured signal from the detection areas, e.g. before step b) in the method for determining 
the amount of a substance, retaining the radiation associated substance before the detection 
areas/particle beds of the disc have been treated, e.g. with antibody-fluorochrom. All search 
areas/detection areas are scanned and each of the corresponding sets of measured data is stored 
in a background radiation image. When using antibody-fluorochrom, the corresponding image is 
called background fluorescence image. The image contains the measured data of the background 
radiation noise signal and other noise signals. In preferred variants the scanning for the 
background radiation image is carried out as close as possibly to the step in which the substance 
causing the radiation is formed and/or retained in the detection microcavities. Although not 
preferred, the manufacturer may provide a background radiation image that for instance may be 
based on the average of many CDs of the same kind. 

The second step, step a'2, is in this first variant the same as step c) in the method 
for determining the amount of a substance, above. The substance is now present in the detection 
microcavities, for instance retained in particle beds that may have been washed, e.g. after the 
particle beds have been treated with e.g. antibody-fluorochrom. This image is the substance raw 
data image, described above. 

If the background radiation images and corresponding substance raw data images 
exist already and are availably stored, there is consequently only to access (e.g. by copying or 
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reading) the images in step a' 1 and step oc'2 instead of generating them again by the scanning 
and detecting processes. 

The two images should be detected under the same conditions, e.g. same laser 
power and same detector sensitivity if the detection principles is based on fluorescence. 

The third step of the first variant, step a' 3, is to correlate the two images in two 
dimensions, that define a correlation image, by use of a two-dimensional correlation. The use of 
two-dimensional correlation is well-known in the art of signal-processing, see reference [1], page 
90-92. Thus, the step has resulted in a correlation image. 

In the fourth step, step oc'4, of the variant, the position of the pixel, or pixels, 
having the largest intensity value in the correlation image is/are located. This could be done by 
using a simple comparison algorithm, wherein the recently largest value is stored and compared 
with the next data, e.g.intensity value, that hasn't been compared with any of the earlier pixel 
data. The positions of the pixels having the largest pixel values, largest intensity values, are 
stored and used in the next step, step a' 5. 

In step a' 5, the background radiation image is moved to the correct position of the 
correlation image by using the identified pixel positions. The parts of the image that overlap is 
subtracted pixel by pixel. 

These five steps have resulted in a new image, the background processed image. 

The substance raw data image comprise the wanted substance radiation data. If not 
any new noise have been introduced by the introduction of the substance giving rise to the 
measured signal into the microcavities (part of step b), the process of noise removing may be 
stopped. 

In the second variant-is a Median value of the background radiation data used. In 
figure 6, the first variant is presented schematically in a flowchart. Following steps may be used: 

accessing or generating a substance raw data image (step a" 1); 

determining a Median value for background radiation data (step of '2); 
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subtracting pixel values of substance raw data image with the median value, 
and store background processed image (step a"3). 

The second variant of reducing background radiation and construction of 
background processed images will now be described in more detail. 

This second variant for removing background radiation, like fluorescence, 
especially for removing electronic and optical noise that exists in the search area outside the 
detection area, is possible to use instead of the five steps described above. This alternative 
variant works very well under the condition that there is no unwanted object of large size outside 
the detection area and the background signal intensity is constant inside the image. This variant 
has the advantage that the disc is scanned only once, which is described hereinafter. 

The first step, step a" 1 in this alternative variant, should start after the radiation 
associated substance has been incorporated into the detection microcavities. In the method for 
determining the amount of a substance above this for instance corresponds to after the detection 
areas have been treated with the antibody-fluorochrom, i.e. in step c. As described above, during 
the scanning of the search areas for the radiation, the pixel values of the substance raw data 
images are detected and determined, and stored as digital pixel data. If the substance raw data 
images exist already and are availably stored, there is consequently only to copy/read the images, 
instead of generating them again by scanning and detecting processes. 

In the second step, step a median value of the pixel values for at least a 

portion of the search area outside the detection area is calculated and used as background 
radiation value, e.g. as a value for the background fluorescence. The median value is an 
estimated value for the background radiation (intensity). 

In step cc"3, the substance raw data is decreased with the median value 

contribution, e.g. by subtraction, pjxel b\ pixel, until all pixels are processed and a background 

\ 

processed image is provided. 

Potentially the step for determining Median value may be replaced by determining 
some other form of an average value and use this value instead of the Median in the subsequent 
steps. 
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Both methods, a' and a'\ provide a background processed image for each scanned 
search area, i.e. an image with reduced background radiation, such as fluorescence, and can be 
used separately or in combination in other steps of the innovative method, with the former being 
preferred at the moment. 



Step P: Reducing peak disturbances. 

The method of this aspect of the invention may comprise one or more steps in 
which peak disturbances originating from dust, chemistry and other artefacts are reduced, e.g. 
reducing peaks from fibres. The disturbance peaks have typically a relatively high amplitude and 
a relatively narrow distribution. However, the image of an ideal detection area is smooth and 
does not contain any disturbance peaks. 



In figure 7, step p for reducing the peak disturbances is presented schematically in 
a flowchart. Following steps may be used: 

- Generating a Laplacian processed image from the earlier stored image, 
typically a background processed image, by using a Laplacian Filter in two 
dimensions (step pi); 

- Calculating a disturbance threshold value from the image created in the 
previous step, typically from a background processed image (step p2); 

Generating binary encoding Laplacian processed image by applying the 
disturbance threshold limit to the image used for calculating the treshold value 
(step p3); 

- Finding the start points in clusters (step P4); 

- Determining a minimum rectangular area surrounding the start points (step p5); 

- Calculating an interpolated value within each peak area and substituting the 
original data of the peaks, and storing the disturbance filtered image (step P6). 
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The process for reducing peak disturbances will now be described in more detail. 

Thus in a first step, step (31, a radiation intensity image, for instance a background 
processed image as discussed above, may be processed by a Laplacian filter in two dimensions. 
This is a standard method earlier known from reference [1]. The result of this step is that peaks 
in a background processed image are amplified. The result is normalized with the maximum 
radiation (intensity) in the background image. This step provides a Laplacian processed image. 

In a second step, step p2, a disturbance threshold limit is calculated by use of the 
standard statistic functions average and standard deviation of the stochastic background noise, 
calculated from the image used in step p 1 . This step enhances the disturbance peaks of the 
*image. 

In the third step, step (33, the Laplacian processed image is binary encoded. In the 
encoding process pixel intensities/values equal to or higher than a threshold limit are set to one 
of the two corresponding binary values, for instance "1" or "positive" and pixel intensities less 
than the disturbance threshold limit are set to the other of the binary values, for instance "(T or 
"negative". The determined disturbance threshold limit determined in step p2 may be used in 
step P3 for encoding. This step creates an image wherein the high intensity pixels of the peaks 
appears as groups, even called clusters. 

In a fourth step, step p4, the Laplacian processed image may then be processed for 
finding the pixel in each peak having the largest value by using the position information from the 
clusters in the binary Laplacian processed image. Each pixel having the largest value in a cluster 
is defined to be a center and start point. The other pixels of the clusters are set low in the binary 
Laplacian processed image. This step results in a cluster image. 

In a fifth step, step p5, the start points defined in the cluster image is used for 
determining a minimum rectangular area by iterative testing using different area sizes. The 
iteration, that works on the Laplacian processed image, is stopped when a minimum rectangular 
area contains only pixels having a negative intensity value or zero value according to the 
definition in previous steps and as a result of the earlier performed Laplacian filtering. Said 
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pixels surround the corresponding start point and defines the position, the area of the 
corresponding disturbance peak and the boarder pixels which is outside the peak area. 

In a sixth step, step 06, the disturbance peaks of the background processed image 
are removed by substituting the original data within each peak area with an interpolated value 
between the boarder pixels of the corresponding peak. The interpolated values are calculated by 
use of an appropriate interpolation function, both a linear and a non-linear function is a possible 
choice. 

The process has now generated a "disturbance filtered image" of the radiation 
intensity, e.g. fluorescence intensity, of the detection area. 



Step x: Determining a global threshold value 

The object of this part of the invention is to make a radiation intensity image, such 
as a disturbance filtered image, binary by determining a global threshold that separates the 
measured substance signals, i.e. the signals from the substance in the detection area, from the 
background signals. For instance, the detection area pixels should be equal to binary high, like 
"1", and the rest of the image equal to binary low, like "0". 

The purpose of the following steps is to determine a valid global threshold value. 
This is done by an iteration that tests different threshold values. A global threshold is valid when 
only a few pixels in a part of the search area that is not connected to the detection area have an 
intensity higher than the pixels in the detection area. 

In figure 8, step % for determining a global threshold value is presented 
schematically in a flowchart. Following steps may be used: 

- Generating a y-summation graph (step % 1 ); 

- Median filtering of the y-summation graph (step x2); 

- Determining centre position of detection area (step ^3); 

- Determining radial directed edges of detection area (step %4); 
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Generating a x-summation graph (step %5); 

Median filtering of the x-summation graph (step %6); 

- Determining an edge curve (correlation) (step x7); 

- Restraining disturbance peaks in the edge curve (step 

5 - Determining the start and end positions of the detection area (step %9); 

Determining a start value of the global threshold value (step %\0); 

- Adjusting the global threshold value to the value giving correct detection area 
width (step 

Generating a binary encoded disturbance filtered image by use of the global 
10 threshold value (step %\2). 



This process will now be described in more detail. 

In a first step, step xl > pixel values are summarized in the radius direction, pixel 
values, line by line, in a radiation intensity image, typically the disturbance filtered image, to get 
1 5 a one-dimensional y-summation graph in the y-direction of the radiation intensity, e.g. 

fluorescence intensity (angular direction on a spinning disc). Each y-coordinate on the graph will 
be the sum of the pixel values in radial direction having the same y-coordinat. 

In a second step, step %2«the y-summation graph is filtered to remove disturbing 
peaks, e.g. with a median filter function. 

20 In a third step, step ^3. the centre of the detection area in y-direction (angular 

direction on a spinning disc) is calculated and located by use of a correlation, and with the 
filtered y-summation graph and the known detection area width as input value to the correlation 
length. This should normally give a graph with one large peak. 



25 



In a fourth step, step the centre of the largest peak is set to the centre of the 
detection area. It's now possible to calculate the positions and the radial extension of the particle 
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bed edges by use of the known data for the centre and detection area/particle bed width. Hereby, 
the radial directed boarder/edge pixels are known. 

In a fifth step, step %5, the radiation intensities, are summarized in the image in y- 
direction (angular direction on a spinning disc) using the limits (boarders) determined in the 
5 previous step to get a x-summation graph (radial summation graph) in x-direction (radius 

direction) of the radiation intensities. The summation of pixel values is performed in the same 
way, but in the orthogonal direction in the same plane, as in step %\ . 

In the next step, step x6, the x-summation graph is filtered with a median filter 
function to remove disturbing peaks. 

10 In the seventh step, step %7, the start of the detection area in the x-direction is 

determined by use of correlation between a step function and the filtered x-summation graph. 
The mathematical process results in an edge curve, see figure 9. 

In step x8, disturbance peaks, close to edge that is sitting opposite to the inlet, are 
restrained by multiplying the edge curve with a ramp function. 

15 The purpose of step y9* is t0 find the start position of the detection area by finding 

the highest positive peak centre in the filtered edge curve (from step %%) and using a minimum 
detection area length to calculate the end position of the detection area. 

An approximate position of the detection area edges, boarder pixel is now known. 

In the tenth step, step xlO, the global threshold value is now determined by use of 
20 optimal threshold technique for image processing (reference [1], starting at page 360). It is 
assumed that the background and particle bed peaks are normal distributed with the same 
variance. Calculate from the disturbance filtered image the estimated background intensity as 
max of (median background intensity on the areas A1+A3 or median intensity background on the 
area A2)). The positions and orientations of the areas Ah A2 and A3 is illustrated in figure 13. 
25 Calculate the estimated substance intensity as median of the detection area. Calculate the start 
value for the iteration step as average of the estimated background intensity and the estimated 
substance intensity. 
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In step %\ I, the global threshold value is iterative adjusted, between a maximum 
intensity value, defined by the estimated substance intensity, and a min value, defined by the 
estimated background intensity, so that the average number of pixels in the y-direction (the 
angular direction) inside the approximate position of the detection area is equal to the detection 
5 area width. 

In step xl2, the calculated global threshold value is used for binary encoding a 
stored image, in particular the stored disturbance filtered image into a binary coded disturbance 
filtered image. 



10 Step 5: Remove binary objects outside the detection area 

The object of the following steps 51 - 52, is to remove any disturbances outside the 
detection area. The following algorithm processes the data of the binary coded disturbance 
filtered image, in which all objects are set to 1, i.e. binary high. Further, it is assumed that there 
is no object to be removed inside the boarders of the approximately calculated detection area 
15 position, and the detection area is the largest area of the image. 

In figure 10, step 5 for removing binary objects outside the detection area is 
presented schematically in a flowchart. 

The first step, step 51 : labelling of all the objects in the binary coded disturbance 
filtered image, wherein the number of connected pixels, i.e. the area, of each object is 
20 determined. Labelling is a known method in image processing. 

The second step, step 52: generating a binary detection area image by removing all 
but the largest objects that are inside the approximate detection area position. 

Step s: Removing unwanted areas of the detection area 

25 Sometimes the object is to calculate an amount for a "normalized'* area, i.e. a 

predefined amount of pixels, inside the detection area. This is done to minimize the variance 
between detection areas on the same microfluidic device. The variance depends among other 
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things on background radiation, such as fluorescence, from the detection area. It is assumed that 
the 3D distribution of an intensity signal inside the detection area is independent of small 
variations in the detection area length. The substance signal deriving from a portion of the 
detection area that may be closests to the inlet of the corresponding detection microcavity. The 
algorithm used is a process for generating a new binary image that only comprises pixels within 
the normalized area ((detection area width) multiplied by (min detection area length)) in the 
binary detection area image starting from the inner part of the CD. 

In figure 11, step e for removing unwanted areas of the detection area is presented 
schematically in a flowchart. 

The first step, step si: generating a normalized area image by setting all binary 
high indicated pixels outside the location of the normalized area ((detection area width) 
multiplied by (min detection area length)) to binary "low" in the binary detection area image. 

The second step, step s2: multiplying each corresponding pixel in the disturbance 
filtered image, generated and stored earlier in the invented process, and the normalized area 
image with each other and perform a calculation, e.g. integration, on the normalized area. 

Step (j>: Apply default detection area in noisy images 

The described method above will not work properly, or not at all, for very noisy 
images hiding the substance signal. The following algorithm works, if all processed images 
originate from the same set of detection areas, e.g. the same CD. The detection areas are then 
well-defined in angular position and quite well in radius direction. 

During the algorithm to locate the position of the detection area within the search 
area, all positions for detection-areas having high confidence are stored. 

In figure 12, step (j) for applying a default detection area in noisy images is 
presented schematically in a flowchart. 

The method for applying default detection area in noisy images may comprise 
following steps: 
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Step §1: Saving automatically all determined detection area positions (step %) 
fulfilling high confidence criterion. 

Step §2: Calculating an average detection area position. 

Step <j>3: Using automatically the average detection area position if one or more 
5 predetermined condition is fulfilled. 

The average position for a detection area will automatically be applied for all 
calculated detection areas that fulfil all following conditions: 

• It's angular start position differ from the average start position more than 
[amount] multiplied by [detection area width], wherein the [amount] is chosen 0-1. 

10 • The average number of binary high pixels, in the normalized area image, 

in angular direction for all radiuses is less than [amount] multiplied by [detection area width], 
wherein the [amount] is chosen 0-1. 

• The difference between angular stop and start positions is greater than 
[amount] multiplied by [detection area width], wherein the [amount] is chosen 1-10. 

15 • The number of binary high pixels, in the normalized area image, in 

angular direction for all radiuses is less than [amount] multiplied by [detection area width] 
multiplied by [min detection area length], wherein the [amount] is chosen 0-1. 

Steps d) and e) in the method for determining the amount of a substance that is 
20 associated with the collected radiation as well as any of the other individual method steps or their 
combinations discussed are typically performed by the appropriate software lor instance included 
in the controller or elsewhere, for instance not in physical association with the innovative 
arrangements described herein. 

The invention also comprises computer program-related aspects for treating radiation 
25 data that have been assembled by steps (a)-(c) of the innovative method. One such aspect is a 
computer program product that 
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(1) comprises program code elements corresponding to a sequence variant of step (d) 
above comprising one or more of the substeps described for this step, and 

(2) when installed on the appropriate hardware is capable of causing the hardware 
(computer) to execute the sequence of substeps on data which have been obtained by 

5 performing step (c) above on any kind of microfluidic device irrespective of being 

spinnable or not. 

Other code elements may be included, for instance corresponding to step (e) in order 
to execute the sequence (d)-(e). Another computer program-related aspect is the computer 
program product stored at a computer program readable means which, when the product is 
10 loaded, makes it possible for a computer to perform the sequence of steps corresponding to the 
code elements in the stored computer program product. A third computer program- related aspect 
is a carrier having at least one of the innovative computer program products thereon. The carrier 
may be a computer memory, a Read-Only Memory or an electrical signal carrier. 

The present invention may be implemented as a computer program product directly loadable into 
15 an internal memory storage of a processing unit within the controller in the detector arrangement, 
comprising the software code means for performing the steps of any claims 1-16. 

Further, the present invention relates to a computer program product stored on a computer usable 
medium, comprising readable program for causing a processing unit in a computer means to 
20 control an execution of the steps of any of the claims 1-16. 

The computer usable medium is a record medium, a computer memory, a Read-Only Memory or 
an electrical carrier signal. 

25 Certain innovative aspects of the invention is defined in more detail in the appending 

claims. Although the present invention and its advantages have been described in detail, it 
should be understood that various changes, substitutions and alterations can be made herein 
without departing from the spirit and scope of the invention as defined by the appended claims. 
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Moreover, the scope of the present application is not intended to be limited to the particular 
embodiments of the process, machine, manufacture, composition of matter, means, methods and 
steps described in the specification. As one of ordinary skill in the art will readily appreciate 
from the disclosure of the present invention, processes, machines, manufacture, compositions of 
matter, means, methods, or steps, presently existing or later to be developed that perform 
substantially the same function or achieve substantially the same result as the corresponding 
embodiments described herein may be utilized according to the present invention. Accordingly, 
the appended claims are intended to include within their scope such processes, machines, 
manufacture, compositions of matter, means, methods, or steps. 
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CLAIMS 

1. Method for reducing noise in an image, typically a substance raw data image, 
comprising substance and noise data which is associated to one detection microcavity 
and one substance is present in the detection microcavity of a microchannel e 

5 structure of a microfluidic device, the substance raw data information obtained by use 

of a detector arrangement capable of scanning and collecting radiation data in a 
search area associated to a detection area, the collected radiation data are stored as 
measure data for the amount of the substance as substance raw data in a set in a data 
storage, said method is characterized in that the method comprises following step: 
10 - reducing background radiation data either by use of a image comprising 

background radiation data and no substance data information, a background radiation 
image, (step a') or, by use of a calculated median value of the background radiation 
data (step a"). 

2. Method according to claim 1, characterized in that the method comprises the step of: 
15 - if available, accessing a background radiation image or, if not available, generating 

a background radiation image (step a'l); 

- if available, accessing corresponding substance raw data image or, if not available, 
generating a substance raw data image (step a '2); 

- correlating the two images and generating a correlation image (step a 5 3); 
20 - locate the positions of the largest pixel values (step of 4); 

- positioning the images correctly, subtracting corresponding pixel values and storing 
the generated background processed image(step a' 5). 

3. Method according to claim 1, characterized in that the method comprises the step of: 

- if available, accessing corresponding substance raw data image or, if not available, 
25 generating a background radiation image substance raw data image (step a" 1); 

- determining a Median value for background radiation (step a' '2); 

- subtracting pixel values of substance raw data image with the median value 
contribution, and store background processed image (step ct"3). 

4. Method according to claim 1, characterized in that the method comprises the step of: 
30 Reducing peak disturbances (Step (3), 
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5. Method according to claim 4, characterized in that the step of reducing peak 
disturbances (Step (3) comprises following steps: 

- Generating a Laplacian processed image from the earlier stored image, typically a 
background processed image, by using a Laplacian filter in two dimensions (step (31); 

5 - Calculating a disturbance threshold value from the earlier stored image, typically a 

background processed image (step (32); 

- Generating binary encoding Laplacian processed image by use of the disturbance 
threshold limit (step (33); 

- Finding the start points in clusters (step (34); 

10 - Determining a minimum rectangular area surrounding the start points (step (35); 

- Calculating an interpolated value within each peak area and substituting the original 
data of the peaks, and storing the disturbance filtered image (step P6). 

6. Method according to claim L characterized in that the method comprises the step of: 
Determining a global threshold value (step x)- 

1 5 7. Method according to claim 6, characterized in that the step of determining a global 

threshold value (step x) comprises following steps: 

- Generating a y-summation graph (step 

- Median filtering of the y-summation graph (step x2); 

- Determining centre position of detection area (step x3); 
20 - Determining radial directed edges of detection area (step 

- Generating a x-summation graph (step x5); 

- Median filtering of the x-summation graph (step x6); 

- Determining an edge curve (correlation) (step x7); 

- Restraining disturbance peaks in the edge curve (step x8); 

25 - Determining the start and end positions of the detection area (step x9); 

- Determining a start value of the global threshold value (step xlO); 

- Adjusting the global threshold value to the value giving correct detection area width 

(stepxll); 

- Generating a binary encoded disturbance filtered image by use of the global 
30 threshold value(step y\2). 
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8. Method according to claim 1, characterized in that the method comprises the step of: 
removing binary objects outside the detection area (step 5). 

9. Method according to claim 8, characterized in that the step of removing binary 
objects outside the detection area (step 5) comprises following steps: 

5 - labelling of all the objects in the binary coded disturbance filtered image, wherein 

the number of connected pixels of each object is determined (step 51); 

- generating a binary detection area image by removing all but the largest objects that 
are inside the approximate detection area position (step 52). 

10. Method according to claim 1, characterized in that the method comprises the step of: 
10 removing unwanted areas of the detection area (step s). 

11. Method according to claim 10, characterized in that the step of removing unwanted 
areas of the detection area (step s) comprises following steps: 

- generating a normalized area image by setting ail binary high indicated pixels 
outside the location of the normalized area ((detection area width) multiplied by (min 

15 detection area length)) to binary "low" in the binary detection area image (step si); 

- multiplying each corresponding pixel in the disturbance filtered image, generated 
and stored earlier in the invented process, and the normalized area image with each 
other and perform the calculation, e.g. integration, on the normalized area (step s2). 

12. Method according to claim L characterized in that the method comprises the step of: 
20 Applying default detection area in noisy images (step 

13. Method according to claim 12, characterized in that the step of applying default 
detection area in noisy images comprises following steps: 

- Saving automatically all detection area positions fulfilling high confidence criterion 

(step cj) 1 ); z 
25 - Calculating an average detection area position (step §2): 

- Using automatically the average detection area position if one or more 
predetermined condition is fulfilled (step <|>3). 

14. A method for determining a measure of at least one substance comprising the steps 
of: 
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(a) providing a micro fluidic device and a detector arrangement, wherein 
said microfluidic device comprises a plurality of microchannel structures, each 
of which has an inlet port and a detection microcavity, and a plurality of 
detection areas, each of which is associated with one of said detection 

5 microcavities; and wherein said detector arrangement is capable of collecting 

radiation from individual subareas of each of said detection areas; 

(b) processing one or more liquid aliquots in at least one of said plurality 
of microchannel structures so that said substance is retained in the detection 
microcavity of each of said at least one of said plurality of microchannel 

10 structures; 

(c) scanning the detection areas associated with the detection 
microcavities that are part of microchannel structures in which step (b) has been 
carried out to obtain radiation from individual subareas (pixels) of each scanned 
detection area, said scanning being performed by the use of said detector 

1 5 arrangement; 

(d) integrating radiation as a function of subarea for each scanned 
detection area to obtain the amount of radiation from each detection area; and 

(e) characterizing for each of the amounts obtained in step (d) a reaction 
variable that has been included in the process protocol used for each 

20 microchannel structure. 

15. The method of claim 14, wherein at least one of the steps in claims 1-13 is carried out 
in step (d). 

16. The method of claim 15, wherein integration only is on a portion of each detection 
area, a normalized ajea. 

25 17. A computer program product directly loadable into an internal memory storage of a • 

processing unit within the computer means, such as a controller of a detector 
arrangement, comprising the software code means for performing the steps of any 
claims 1-13 and/or 14-16. 
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18. The computer program product stored on a computer usable medium, comprising 
readable program for causing a processing unit in a computer means, such as a 
controller of a detector arrangement to control an execution of the steps of any of 
claims 1-13 and/or 14-16. 

5 19. The computer program product according to claim 18, characterised in that the 

computer usable medium is a record medium. 

20. The computer program product according to claim 18, characterised in that the 
computer usable medium is a computer memory. 

21. The computer program product according to claim 18, characterised in that the 
10 computer usable medium is a Read-Only Memory. 

22. The computer program product according to claim 18, characterised in the computer 
usable medium is an electrical carrier signal. 

23. Arrangement for reducing noise in an image, preferably a substance raw data image, 
comprising substance and noise data which is associated to one detection microcavity 

15 and one substance retained in the detection microcavity of a micro fluidic structure of 

a microfluidic device, the substance raw data information obtained by use of a 
detector arrangement capable of scanning and collecting radiation data in a search 
area associated to a detection area, which is associated to the detection microcavity, 
the collected radiation data are stored as measure data for the amount of the substance 

20 as substance raw data in a set in a data storage, the arrangement also comprises a 

processing unit for controlling the arrangement, characterized in that the method 
comprises following step: 

- means for reducing background radiation data either by use of an image 
comprising background radiation data and no substance data information, typically a 

25 background radiation image, (step a') or, by use of a calculated median value of the 

background radiation data (step a"). 

24. Arrangement according to claim 23, characterized in that the arrangement comprises: 
- means for accessing a background radiation image and means for generating a 
background radiation image (step of 1); 
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- means for accessing corresponding substance raw data image and means for 
generating a substance raw data image (step a' 2); 

- means for correlating the two images and generating a correlation image (step a' 3); 

- means for locate the positions of the largest pixel values (step oc'4); 

5 - means for positioning the images correctly, means for subtracting corresponding 

pixel values and means for storing the generated background processed image(step 
ot'5). 

25. Arrangement according to claim 23, characterized in that the arrangement comprises: 

- means for accessing corresponding substance raw data image and means for 
10 generating a background radiation image substance raw data image (step a' " 1 ); 

- means for determining a Median value for background radiation data (step a'*2); 

- means for subtracting pixel values of substance raw data image with the median 
value contribution, and means for storing the background processed image (step 
a' ? 3). 

15 26. Arrangement according to claim 23, characterized in that the arrangement comprises 

means for reducing peak disturbances (Step p). 

27. Arrangement according to claim 26, characterized in that the means for reducing peak 
disturbances (Step P) comprises: 

- means for generating a Laplacian processed image from the earlier stored 

20 background processed image by using a Laplacian filter in two dimensions (step pi); 

- means for calculating a disturbance threshold value from the background processed 
image (step P2); 

- means for generating binary encoding Laplacian processed image by use of the 
disturbance threshold limit (step (J3); 

25 - means for finding^the start points in clusters (step p4); 

- means for determining a minimum rectangular area surrounding the start points 
(step P5); 

- means for calculating an interpolated value within each peak area and substituting 
the original data of the peaks, and storing the disturbance filtered image (step p6). 
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28. Arrangement according to claim 23, characterized in that the arrangement comprises 
means for determining a global threshold value (step x). 

29. Arrangement according to claim 28, characterized in that means for determining a 
global threshold value (step x) comprises: 

5 - means for generating a y-summation graph (step xl); 

- means for Median filtering of the y-summation graph (step %2)\ 

- means for determining centre position of detection area (step x3); 

- means for determining radial directed edges of detection area (step x4); 

- means for generating a x-summation graph (step x5); 

10 - means for median filtering of the x-summation graph (step y6)\ 

- means for determining an edge curve (correlation) (step yl)\ 

- means for restraining disturbance peaks in the edge curve (step x8); 

- means for determining the start and end positions of the detection area (step x9); 

- means for determining a start value of the global threshold value (step xlO); 

15 - means for adjusting the global threshold value to the value giving correct detection 

area width (step xl 1); 

- means for generating a binary encoded disturbance filtered image by use of the 
global threshold value (step x!2). 

30. Arrangement according to claim 23, characterized in that the arrangement comprises 
20 means for removing binary objects outside the detection area (step 8). 

3 1 . Arrangement according to claim 30, characterized in that the means for removing 
binary objects outside the detection area (step 5) comprises: 

-means for labelling of all the objects in the binary coded disturbance filtered image, 
wherein the number of connected pixels of each object is determined (step 81 ): 
25 - means for generaSng a binary detection area image by removing all but the largest 

objects that are inside the approximate detection area position (step 82). 



32. Arrangement according to claim 23, characterized in that the arrangement comprises 
means for removing unwanted areas of the detection area (step s). 
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33. Arrangement according to claim 32. characterized in that means for removing 
unwanted areas of the detection area (step s) comprises: 

- means for generating a normalized area image by setting all binary high indicated 
pixels outside the location of the normalized area ((detection area width) multiplied 
by (min detection area length)) to binary "low" in the binary detection area image 
(step si); 

- means for multiplying each corresponding pixel in the disturbance filtered image, 
generated and stored earlier in the invented process, and the normalized area image 
with each other and perform the calculation, e.g. integration, on the normalized area 
(step s2). 

34. Arrangement according to claim 23. characterized in that the arrangement comprises 
means for applying default detection area in noisy images (step (j)). 

35. Arrangement according to claim 34, characterized in that the means for applying 
default detection area in noisy images comprises: 

- means for having automatically all detection area positions fulfilling high 
confidence criterion (step (|)1); 

- means for calculating an average detection area position (step (|)2); 

- means for using automatically the average detection area position if one or more 
predetermined condition is fulfilled (step <j>3). 

36. An arrangement for determining a measure of at least one substance, the arrangement 
comprises: 

(a) a microfluidic device and a detector arrangement, wherein said 
microfluidic device comprises a plurality of microchannel structures, each of 
which has an inlet port and a detection microcavity, and a plurality of detection 
areas, each ofi\vhich is associated with one of said detection microcavities: and 
wherein said detector arrangement is capable of collecting radiation from 
individual subareas of each of said detection areas: 

(b) means for integrating radiation as a function of subarea for each 
scanned detection area to obtain the amount of radiation from each detection 
area: and 
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(c) means for characterizing for each of the amounts obtained in step (d) 
a reaction variable that has been included in the process protocol used for each 
microchannel structure. 

37. The arrangement according to claim 36, at least one of the method steps a - § in 
5 claims 1-13 is carried out in step (d). 

38. The arrangement according to claim 37, characterized in that the subarea in step 
d) is a normalized area. 
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